BI-TP 2008/09 



Baryon number, strangeness and electric charge 
fluctuations at zero and non-zero chemical potential 

C. Schmidt (for the RBC-Bielefeld Collaboration) 

Universitat Bielefeld. Fakultat fur Physik 
^ ! Postfach 100131, D-335G1 Bielefeld, Germany 



O 

(N 



(N 



> 

cn 
(N 
O 



E-mail: schmidtfiphysik . uni-bielef eld . de 



' Abstract. We present results on baryon number, strangeness and electric charge 

fluctuations in QCD at non-zero density and temperature obtained from lattice 
calculations with almost physical quark masses. At vanishing chemical potential, 
i.e. under conditions almost realized at RHIC and the LHC, quartic fluctuations 
of net baryon number and strangeness are large in a narrow temperature interval 
' characterizing the transition region from the low to the high temperature phase. Our 

results are based on Taylor expansions in light and strange quark chemical potentials, 
Oh| i.e. we rigorously compute corrections to bulk thermodynamic quantities at non 

vanishing chemical potential, by performing a Taylor expansion in fi/T. We find 
non-monotonic behavior for the radius of convergence of this series, which could be a 
hint for a critical end-point in the (T, /i)-plane. 



1. Introduction 



O ' At non-zero chemical potential, lattice QCD is harmed by the "sign-problem", which 
00 ' 

^ I makes direct lattice calculations at non-zero density practically impossible by all 
known lattice methods. Some methods, however, allow to extract information on 
I the dependence of thermodynamic quantities on a small chemical potential parameter, 
H ' based on lattice calculations performed at zero or imaginary chemical potential. For an 
overview see, e.g. [1]. 

We will report here preliminary results on the Taylor expansion method with almost 
realistic quark masses, i.e. a realistic strange quark mass and a light quark mass which 
is about a factor of 2 heavier than the physical quark mass. The analysis is based on 
data generated in the context of the recent equation of state calculation of the RBC- 
Bielefeld Collaboration [2j . The logarithm of the grand canonical partition function or 
equivalently the pressure (p) can be Taylor expanded in /iu/T, /i^/T and yU^/T where 
/^n,d,s are the up-, down- and strange-quark chemical potentials, respectively, and T is 
the temperature. We define 

(1) 

The expansion coefficients c^f^iT) are computed on the lattice at zero chemical 
potential, using stochastic estimators. For details see [3]. 
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Figure 1. Quadratic (left) and quartic fluctuation (right) of baryon number {B), 
electric charge (Q) and strangeness (S), normalized by their corresponding Stefan- 
Boltzmann value. The results are from iVt = 4 lattices, and in the case of strangeness 
fluctuations, open symbols show results from Nt ~ 6 lattices. Thin solid lines are from 
Ferrenberg-Swendsen analyses. 



2. Quadratic and quartic fluctuations 

The QCD partition function is naturally formulated in terms of quark-fields, and thus 
quark chemical potentials. However, alternatively to the quark chemical potentials one 
can introduce chemical potentials for the conserved quantities baryon number B, electric 
charge Q and strangeness S {fiB,Q,s), which are related to fiu,d,s via 

12 11 11 

By means of these relations the coefficients cf'^^^^ of the pressure expansion in terms of 
fJ'B,Q,s are easily obtained, in analogy to Eq. [1] 

|i-E4rm(^)*(W(f)'- (3) 

The quadratic (Gaussian) fluctuations (xf ''^'^) of B, Q and S respectively are related 
by the fluctuation dissipation theorem to the second derivatives of the partition function 
with respect to the corresponding chemical potentials fiB,s,Qy whereas the quartic (non- 
Gaussian) fluctuations (xf ''^'^) are defined by the fourth derivatives. In terms of the 
expansion coefficients we find at zero chemical potential 
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In Fig.[T]we show these fiuctuations as a function of the temperature. The quadratic 
fiuctuations rapidly increase at a transiton temperature T = and quickly approach 
1 for T > Tc. Note, that all fluctuations have been normalized by their corresponding 
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Figure 2. Ratio of quartic and quadratic fluctuations for baryon number B (left) and 
electric charge Q (right). Shown are the preliminary results for 2+1-flavor QCD with 
TTiTr w 220 MeV and the previous results for 2-flavor with ttItt 770 MeV. Thin solid 
line (left) is a Ferrenberg-Swendsen analysis. 



Stefan-Boltzmann value. The quartic fluctuations show a peak at T = T^, which is most 
pronounced for the baryon number fluctuations and least pronounced for the strangeness 
fluctuations. In fact, by using an appropriate scaling Ansatz for the free energy, one 
can show that X4 will develop a cusp in the chiral limit, where the transition becomes 
2nd order. For the strangeness fluctuations we show results from two different lattice 
spacings {Nt = 4, 6) and find a relatively small cut-off dependence, which is of similar 
magnitude as we have found for the pressure itself [2]. 

In Fig. [2] we show the ratios of X2 and X4; for fluctuations of B and Q. Ratios are 
well suited quantities to compare with experiment, since many systematic errors as well 
as the volume dependence are eliminated by taking the ratio. In Fig. [2], we also compare 
our preliminary (2-|-l)-flavor data with a pion mass of m.^ ~ 220 MeV to the previous 
2-flavor lattice results with ^ 770 MeV [1], as well as the resonance gas predictions 
(for T < Tc). Besides a smaller Tc, we find that the fluctuations now significantly rise 
above the resonance gas level, which was previously not observed. Note also, that this 
quantity directly gives access to the relevant degrees of freedom, thus already above 
1.5Tc the quantum numbers of the medium are those of a quasi free gas of quarks and 
gluons [4]. 

3. The radius of convergence and the phase diagram 

Estimating the radius of convergence of the series (Eq. [3]) is a method to locate the 
critical end-point (CEP) in the (T, /iB)-phase diagram. The radius of convergence (p) 
is defined as 

p = lim pn with pn = Jc^;of/Cn+2,o.o (4) 

For T > TcEP, the estimators for the radius of convergence should be large and 
eventually will diverge. For T < Tcep the radius of convergence will be limited by 
the phase transition line. In Fig. [3] we plot the first approximation for the radius of 
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Figure 3. The approximation p2 of the radius of convergence of the Taylor series of 
the pressure with respect to the baryon chemical potential, as explained in the text. 
Shown are the preliminary results for 2+1-flavor with ttIt^ w 220 MeV and the previous 
results for 2-flavor with 111^^ ?» 770 MeV, as well as the resonance gas values for p2 , p4 
and earlier calculations of the CEP [5, 6]. Thin dashed and solid lines indicate a 
suggested phase diagram and are only meant to guide the eye. 



convergence p2 in the (T, /x^) -plane, together with previous results for p2 as well as 
previous estimates for the CEP from [5], |6] • Also shown are the resonance gas values 
of p2 and p4, which seems to be aproached by the lattice data for T<0.85Tc. For 
an undoubtful determination of the critical point higher approximations for p will be 
required, however the non monotonic behaviour which is now seen for m^r ~ 220 MeV 
and has not been seen for ^ 770 MeV might be a first sign of the critical region of 
the CEP. In addition, thin dashed and solid lines indicate a suggested phase diagram 
and are only meant to guide the eye. 
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